Preprint typeset using I^TjrjX style cmulatcapj v. 5/2/11 
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ABSTRACT 

We present an analysis of all XMM-Newton and Suzaku X-ray spectra of the nearby luminous Seyfert 
galaxy Fairall 9. Confirming previous analyses, we find robust evidence for a broad iron line associated 
with X-ray reflection from the innermost accretion disk. By fitting a spectral model that includes 
a relativistically ionized reflection component, we examine the constraints on the inclination of the 
inner accretion disk and the black hole spin, and the complications introduced by the presence of a 
photoionized emission line system. Employing multi-epoch fitting, we attempt to obtain robust and 
concordant measures of the accretion disk parameters. We also clearly see a soft X-ray excess in 
Fairall 9. During certain epochs, the soft excess can be described with same disk reflection compo- 
nent that produces the iron line. However, there are epochs where an additional soft component is 
required. This can be attributed to cither an additional highly-ionized, strongly blurred disk reflection 
component, or a new X-ray continuum component. 

Subject headings: galaxies: individual(Fairall 9) - X-rays: galaxies - galaxies: nuclei - galaxies: Seyfert 
-black hole physics 



1. BLACK HOLE SPIN CONSTRAINTS FROM AGN 

X-ray spectroscopy provides one of the best tools for 
studying the properties of astrophysical black holes. The 
hard X-ray tail in active galactic nuclei (AGN) and 
Galactic Black Hole Binaries (GBHBs) is produced by 
the inverse Compton scattering of soft thermal pho- 
tons from the accretion disk by hot/energetic electrons 
in the disk corona. Tremendous insight can be gained 
by studying the reflection/reproces sing of these coro- 
nal X-rays from the accretion disk (Fabian et al.1 Il989l : 
iRevnolds fe Nowak|[200l |Mijler| [2007|) . Identifying and 
characterizing the disk reflection features in the X-ray 
spectrum enables us to study the physical properties 
of the inner accretion disk as well as the spin of the 
black hole itself. This is possible because the reflection 
spectrum is affected by strong relativistic effects, specifi- 
cally Doppler shift and gravitational redshift which give 
all emission line features a characteristic skewed profile. 
These line profiles, especially that of the strong and iso- 
lated iron-Ka fluorescent line, can reveal the location of 
the innermost stable circular orbit (ISCO) and hence the 
black hole spin. 

Along with its mass, spin is the most fundamental pa- 
rameter describing a black hole and has a significant in- 
fluence on physics close to the black hole. Importantly, 
in an AGN or GBHB, black hole spin is the only major 
source of energy in the system other than accretion and is 
often cited as the likely pow er source for relativistic jets 
(Blandfor d fe ZnaiekHl977] ). Moreover, the distribution 
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of spin across a population of black holes provides a new 
way of probing t he formation an d growth of supermas- 
sive black holes (fV r olonteril l20I0'). However, in order to 
examine these astrophysical consequences of black hole 
spin, it is crucial that we have robust ways of measuring 
spin. 

In the AGN realm, there are now a number of objects 
for which the supermassive black hole (SMBH) spin has 
been meas ured using X-ray reflectio n; for a summary, see 
Table 2 in I Brenneman et all (|2011ft . Many of these pre- 
vious studies have considered objects that show strong 
warm absorbers in their soft X-ray spectra. This has 
led some authors to suggest that the detection of rel- 
ativistically broadened iron lin es is an artifact of incor- 
rectly modeling the absorption ([Miller. Turner fe Reeved 
2008). While argu ments exist against t hese absorption 
dominated models (Reynolds et al. 2009), it is clearly in- 
teresting to study relativistic disk reflection in objects 
without any obvious absorption. Therefore, in this pa- 
per, we study one of the brightest "bare" Seyfert galaxies 
to show relativistic disk reflection but no warm absorber 
- Fairall 9. 



Fairall 9 is a well-studied Seyfert I galaxy (z = 0.047) 
that has never show n any imprints of a warm absorber in 



its X-ray spectrum I 


Revnoldslll997llGondoin et al.H200I: 


Schmoll et al.ll2009k 


EmmanouloDOulos et al.ll201ll). The 



Galactic column in the direction of Fairall is also mod- 



est (3.15 x 10 20 cm- 2 ; lKalberla et al.|[2005h . giving us an 
unobscured view of its soft X-ray continuum. Fairall 9 
shows limit ed X-ray variability on bot h long and short 
timescales (|Markowitz fe Ed clson 200l|) with the excep- 
tion of occasional sudden dips that last one or two 
weeks (jLohfink et all 12012ft . Studies with ASCA and 
early XMM-Newton observation revealed a prominent 
narro w, neutra l iron Kg line and a possibl y broad iron 
line ([Reynolds! Il997t iGondoin etaLl 1200 lft. An anal- 
ysis o f the 2007 Suzaku observation bv iSchmoll et al.l 
(|2009ft confirms the presence of a moderately weak broad 
iron line and estimates the dimensionless spin parame- 
ter of the black hole to be a — 0.6 ± 0.07 under the 
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assumption that the emissivity profile of the X-ray re- 
flection is an unbroken power law with an emissivity 
index constrained to be less than 5. More recently, 
lEmmanoulopoulos et al.1 ()2011h analyzed a moderately 
deep 2010 XMM-Newton pointing (130 ks), and found 
a weak spin constraint (a = 0.39^o' 3 o) that spans the 
range from slowly spinning to rather rapid prograde spin. 
Worryingly, the disk inclination of i = 64±gdeg de- 
rived by lEmmanoulopoulos et al.l (|2011f) is strongly in- 
co nsistent with t h e valu e of i = 44 ± 1 deg measured 
bv iSchmoll et al.l (|2009t l. Most recently, iPatrick et all 
have examined a deep 2010 Suzaku pointing ob- 
tained as part of the Suzaku AGN Spin Survey Key 
Project (PI Reynolds) . They found the spin to be uncon- 
strained when the spectrum is fitted by a blurred ionized 
disk spectrum, and determined yet a third value for the 
disk inclination, i — 33^ deg. 

In this Paper, we aim to bring clarity to the question 
of the black hole spin and disk inclination in Fairall 9 by 
considering all available XMM-Newton & Suzaku point- 
ings, including the long XMM-Newton pointing of 2010 
and the Suzaku Key Project pointing of 2010. To ob- 
tain robust measurements of the accretion disk parame- 
ters and black hole spin it is crucial to understand the 
continuum and its variation to the highest precision pos- 
sible. We highlight the importance of multi-epoch fitting 
techniques in breaking a modeling degeneracy resulting 
from the presence of narrow ionized iron lines, as well as 
uncovering evidence for a soft excess above and beyond 
that resulting from soft X-ray reflection from the a uni- 
formly ionized accretion disk. Using this methodology, 
we derive concordant values for black hole spin and disk 
inclination. 

The outline of this paper is as follows. First, we de- 
scribe the datasets used in this work and briefly discuss 
data reduction techniques ($2]). After a preliminary in- 
vestigation of spectral fits to the individual data epochs, 
we perform a multi-epoch analysis in order to produce 
robust and concordant constraints on on the spin of the 
black hole and the inclination of the disk ([J3]). We end 
with a brief discussion of the lessons learned concerning 
methodologies for extracting black hole spin from X-ray 
spectra (Q. 

2. OBSERVATIONS & DATA REDUCTION 
2.1. Data Reduction 
2.1.1. Suzaku 

A summary of the observations discussed in this paper 
is presented in Table [T] All Suzaku data reduction was 
performed within Heasoft v6.11. The version 20110608 of 
the relevant Suzaku CALDB files were used. The Suzaku 
data were reduced following the standard procedure de- 
scribed in the " Suzaku ABC Guid«0' . 

For the Suzaku X-ray Imaging Spectrometer (XIS; 
iKoyama et alJl2~007t ) , all data were taken in the full win- 
dow mode. The XIS data files were first reprocessed 
as recommended by the ABC Guide. For all operat- 
ing XIS detectors (XIS0, XIS1, XIS3), the source ex- 
traction region was chosen to be a circle of radius of 
~ 4 arcmin and the background was taken from a region 



on the same chip. Individual spectra and response files 
for each detector and data mode were created using the 
tasks xisrmf gen and xissimarf gen. Finally, the spectra 
and response files of all data modes for a given detector 
were summed, weighting by exposure time, to yield one 
spectrum for each XIS. Light curves were also created 
from the reprocessed event files using x select. 

Th e Hard X-ray Detector (HXD; iTakahashi et al.1 
|2007|) data were first reprocessed as recommended by 
the ABC Guide, and then filtered using xselect with 
the standard criteria. Only the PIN data is considered 
in this work as GSO is not suitable for sources with 
the count rate as low as Fairall 9. For the PIN, the 
appropriate "tuned" background files were downloaded 
from the High Energy Astrophysics Science Archive Re- 
search Center (HEASARC). We then extracted source 
and background spectra from the "cleaned" event files 
using hxdpinxbpi. 

In addition to the reduction steps outlined above, the 
May-2010 Fairall 9 Suzaku data had to be treated with 
special caution as it is affected by a variety of calibra- 
tion/pointing problems. As one can see in Fig. [TJ the 
pointing was unstable and the standard attitude solution 
was inaccurate during this observation, with the satellite 
wobbling back and forth between two positions separated 
by 1 arcmin. Although it is known that Suzaku is suffer- 
ing from an attitude control problem and therefore a mild 
form of this is already accounted for in the calibration 
files delivered with the data, it appeared to be too strong 
to be properly removed with the standard correction 
(Suzaku memo 2010-05). The best results are achieved 
by creating a new, corrected attitude file using the tool 
aeattcor . si written by Mike N owak; a detailed de - 
scription of this tool can be found in lNowak et al.1 (|2011l ). 

Once corrected for, the wobble leaves three remain- 
ing imprints on the spectra. Firstly, since the wobble 
puts some of the source photons into a region of the im- 
age plane where the high-energy vignetting function is 
more poorly known, there are some deviations between 
the different XIS spectra above 8-9 keV. We choose to 
use XIS data up to lOkeV but must be cautious about 
interpreting spectral residuals at the highest XIS ener- 
gies. Secondly, as is well known, the XIS spectra are 
contaminated by absorption from a hydrocarbon layer 
residing on the optical blocking filter. This contamina- 
tion is monitored by the Suzaku team and corrected for 
in the construction of the effective area file. However, 
the contamination has spatial dependence and, for XIS0, 
it appears that the wobble during this observation has 
placed substantial numbers of source photons on parts 
of the image plane in which the contamination is poorly 
known. As a result, the XIS0 spectrum deviates from 
the other XIS spectra below 2 keV. This deviation is sig- 
nificant (reaching values as high as 50 %) and forces us 
to ignore XIS0 data below 2.3 keV. Lastly, due to the vi- 
gnetting issue, there is no reason to suspect that the stan- 
dard XIS/PIN cross normalization factor of 1.18 (nomi- 
nal HXD aimpoint; Suzaku memo 2008-06) is appropri- 
ate for this dataset — the suppression of XIS flux by the 
additional vignetting would increase this cross normal- 
ization. We must allow this cross normalization to be a 
free parameter and, indeed, our best fit models suggest 



7 http:/ /heasarc. nasa.gov/docs/suzaku/analysis/abc/abc. html 
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Black hole spin and soft excess of Fairall 9 



3 



Table 1 

Overview of observations and exposures. Superscript 3 denotes that XIS0/XIS1/XIS3 were used. All Suzaku data were obtained after the 

failure of XIS2. 



Observatory 


Instrument 


Date 


ObsID 


Exposure 
[kscc] 


Suzaku (B) 


XIS a 


2010/05/19 


705063010 


191 




PIN 






162 


Suzaku (A) 


XIS a 


2007/06/07 


702043010 


139 




PIN 






127 


XMM-Newton (B) 


EPIC-pn 


2009/12/09 


0605800401 


91 


XMM-Newton (A) 


EPIC-pn 


2000/07/05 


0101040201 


16 




Figure 1. Image for XIS3 from uncorrected Suzaku pointing for 
Fairall 9. 



values of 1.3. 



2.1.2. XMM-Newton 



The two XMM-Newton datasets were reduced using 
the XMM-Newton Science Analysis System (SAS) ver- 
sion 10.0.0 and the CCF calibration file XMM-CCF- 
REL-270. We only consider EPIC-pn data for the anal- 
ysis performed in this paper, as MOS does not have 
enough signal-to-noise in the iron K region to contribute 
significantly to the constraints. The EPIC-pn data were 
reprocessed and screened using the tool epchain. The 
spectra for the EPIC cameras were then extracted using 
the tool evselect, selecting the default grade pattern. 
Source spectra are taken from a circular region and the 
background spectra from a region on the same chip. The 
effective area files were generated using the SAS task 
arfgen, and the redistribution matrices were produced 
using the task rmf gen. 

3. RESULTS 

3.1. Initial Data Exploration and Model Construction 

Unless explicitly stated otherwise, all XIS/EPIC-pn 
spectra are binned to a signal-to-noise ratio of 10. The 
considered energy ranges for fitting are 0.7-1.5 keV and 
2.3-10 keV for Suzaku-XIS, with the energies around 
2 keV being excluded because of known calibration issues 
around the mirror (gold) and detector (silicon) edges. 
The only exception to this is the XIS0 spectrum for the 
newest Fairall 9 Suzaku pointing; as already mentioned 
in §2.1.1l this spectrum suffers from uncorrected contam- 
ination at soft energies and, therefore, we ignore the re- 
gion of this spectrum below 2.3 keV. EPIC-pn spectra 
are used in the 0.5-10 keV band. For all Suzaku point- 
ings, PIN spectra were binned to S/N of 5 and are used 
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Figure 2. Ratio to a simple power law for Fairall 9 for Suzaku 
A [Top Panel], Suzaku B [Upper Middle Panel], XMM A [Lower 
Middle Panel], XMM B [Bottom Panel]. 

in the 16-35 keV band. As our spectral models include 
a convolution with a relativistic transfer function, which 
requires an evaluation of the underlying model outside 
of the energy range covered by the data, we extended 
all energy grids to energies far beyond the upper energy 
limit given by the highest data bin considered in fitting. 
The observations considered in this work span a duration 
of 10 years (Table Q]), and hence give an idea of the long- 
term variability of the object. None of these observations 
show significant short-term (intra-observation) variabil- 
ity. Therefore, our spectral analysis considers only point- 
ing averaged spectra. 

As an initial exploration of these data, we fit a simple 
absorbed power-law to the 2.3-4.5 keV part of the spec- 
trum; Figure [5] shows the ratio of the data in the 0.5- 
lOkeV band to this simple model for all four datasets. 
All datasets show a prominent narrow, fluorescent iron 
Ka line (6.4 keV rest frame) as well as the blend of the 
K/3/Fe XXVI lines (7.07keV/6.96keV rest frame). The 
spectrum shows clear concavity, revealing the presence 
of a soft excess and, possibly, a Compton reflection. The 
two narrow fluorescent iron lines together with the re- 
flection hump are clear evidence for cold reflection in the 
source, either from the torus or from the outer parts of 
the accretion disk. None of the spectra show any signs 
of warm absorption, as already discussed in ijl] Despite 
these similarities, the spectral shape varies significantly 
throughout the years, with the XMM-Newton pointings 
(2000, 2009) catching the source in a lower flux state. 
This variability is even more apparent when plotting the 
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Figure 3. Unfolded spectra (Suzaku A [black], Suzaku B [red], 
XMM A [purple], XMM B [blue]) of Fairall 9 showing the spectral 
variability of the source. The spectra were rebinncd for plotting. 

unfoldecfl spectra (Fig. [3]) . As the signal-to- noise in the 
short (2000) XMM-Newton pointing is significantly lower 
than the other datasets we do not consider this pointing 
for the remaining part of this work. 

Building upon previous wo rks (jSchmoll et al.l I2009t 
lEmmanoulopoulos et al.l 1 2 llh . we construct a multi- 
component spectral model to describe our spectra. The 
underlying primary continuum is described by a power- 
law with photon in dex T. We use the model pexmon 
(|Nandra et al.l I2007T) to describe the cold reflection of 
this powerlaw from distant material. In particular, 
pexmon self consistently models the strength of the nar- 
row Fe Kq, Fe K/3, and Ni Ka lines, the Compton shoul- 
der of the Fe Ka line, and the Compton reflection con- 
tinuum. The strength of the distant reflection is char- 
acterized by the usual reflection fraction R normalized 
such that R = 1 corresponds to a reflector that subtends 
half of the sky as seen from the X-ray source. As some 
parameters for pexmon cannot be determined from our 
fits we fix them to certain values; we fix the inclination 
to 60 degrees (most probable value for an isotropic distri- 
bution) and abundances to solar values. The high-energy 
cutoff of the continuum powerlaw is fixed to 300 keV. 

To describe the ionized reflection associated with the 
inner accretion disk, we used a modified version of 
ref lionx developed bv lRoss fe~ Fabian (2005). Starting 
from the publicly available version, we redefined the nor- 
malization parameter (norm — > norm/£) such as to sta- 
tistically decouple the true flux normalization from the 
ionization parameter £. This definition results in a more 
rapid and robust convergence of the spectral fit. For the 
ionized reflection model, the iron abundance and ioniza- 
tion parameter are allowed to vary freely. The photon in- 
dex of the irradiating continuum, however, is tied to the 
index for the continuum underlying the cold reflection. 
The ionized reflection component is then relativisticall y 
blurred using the model relconv (iDauser et al.l l2Q10h : 
this naturally gives rise to a relativistic iron line, blurred 
Compton reflection hump and a soft excess. The radial 

9 The unfolding is performed by assuming a diagonal response 
matrix not making any assumptions about the description of the 
spectra. 



dependence of the emissivity of the reflection component 
is assumed to have a broken power-law form, breaking 
from r~ qi to r -92 at a radius i?brcak- The inner edge of 
the X-ray reflection regim e is taken to be at the ISCO 
([Reynolds Sz Fabian 200§), and the outer edge was fixed 
to 400 R g . Provided that q2 > 2, the relativistic blurring 
kernel is only weakly dependent upon this outer radius. 
The accretion disk inclination i, and the black hole spin 
a were left as free parameters. In additio n, p hotoelectric 
Galactic absor ption is modeled with TBnetf^ a newer ver- 
sion of TBabs ([Wilms. Allen fe McCrayll2000t l with cross 
sections set to vern and abundances set to wilm. 

A careful examination of the residuals shown in Fig. [3J 
especially for observation Suzaku-B, reveals a line like 
feature between cold-FeKa and cold-FeK/3 lines. While 
some or all of this emission may be associated with the 
blue peak of the broad iron line, another possibility is 
that we are seeing FeKa line emission from FeXXV 
(which produces a complex at 6.7 keV). It therefore seems 
possible that the spectral model discussed above is miss- 
ing an emission line system producing Fe XXV and pos- 
sibly FeXXVI line emission. Hence we include an addi- 
tional photoionized emission component (Fig. [4]); phys- 
ically, this may be an accretion disk wind (warm ab- 
sorber) which is out of our line of sight. The pho- 
toemission is parametrized by ionization parameter and 
norm. It is modeled using a table model calculated using 
xstar2xspec from the XSTAR model photemi J^l. assum- 
ing irradiation by a power law L £ oc e a with a=-l and 
a number density of 10 10 cm -3 . At the ionization stages 
relevant for this work (log£ > 3.6) there is no dependence 
of the spectral shape on column density for the model, 
its value was therefore kept fixed at 10 22 cm" 3 . We note 
that the inclusion of this emission line component this 
does not preclude the possibility that the residuals at 
6.7 keV are from the blue peak of the iron line, but it 
does allow the model to explore any degeneracy created 
by the superposition of the broad iron line and the nar- 
row ionized line emission. 

When applying this spectral model to the data, a cross 
calibration constant was introduced for Suzaku between 
XIS and PIN (and in case of the 2010 Suzaku pointing, 
also between the XISs). The constant was fixed to the 
values given in the Suzaku ABC Guide, 1.16 for XIS aim- 
point and 1.18 for HXD aimpoint. An exception to this 
is observation Suzaku-B; as mentioned above, due to the 
pointing issues, there is no expectation that the fiducial 
XIS/PIN cross normalization should be appropriate and 
hence it was left as a free parameter. 

We begin with this model as our base model and, in 
£13.21 we fit this model to each dataset individually. After 
noting problematic issues associated with these fits, we 
perform simultaneous fits to data from different obser- 
vations of Fairall 9, under the assumption that the most 
fundamental parameters of the system (spin, accretion 
disk inclination and iron abundance) do not change — 
this so-called multi-epoch fitting is presented in £|3.3l and, 
as we shall see, yields insights into the spectral properties 
of this AGN. 

Throughout this work, the spectral analysis were per- 
formed with the Interactive Spectral Interpretation Sys- 

10 http: / /pulsar. sternwarte.uni-erlangen.de/wilms/research/tbabs/ 

11 http:/ /heasarc. nasa.gov/docs/software/xstar/xstar. html 
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tenO (ISIS Version 1 .6.0-7; [Houck fc Denicolall20q(fl an d 
the newest XSPEC 12.0 models are used (|Arnaudlll996j) . 
All uncertainties are quoted at the 90 % confidence level 
for one interesting parameter (A^ 2 = 2.7). 
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Figure 4. Unfolded spectra (rebinned for plotting), assuming a 
diagonal response, and model [black] of the iron line region for 
observation Su zaku -A fitted with the model including photoionizcd 
emission (see £13,111 . The plot shows the energies in the observed 
frame. 



3.2. Independent fits to individual datasets 

We begin by conducting independent fits of our base 
spectral model to the three remaining observations 
under consideration, XMM-Newton-B, Suzaku A and 
Suzaku B. Our base model provides a good description 
of data well for three datasets, with reduced-^ 2 values 
between 1.00-1.04 (Table [5]). A sample spectrum with 
residuals is shown in Fig. [5] illustrating the good quality 
of the fit. 

We argue below that these fits are incomplete and 
hence that the inferred parameters are thus not neces- 
sarily correct. However, it is instructive to compare our 
results from these fits (Table [5]) with previously pub- 
lished works. For Suzaku A, our derived param eters are 
similar to those found by iSchmoll et al.l (|2009l ) with the 
exception of black hole spin , disk inclination, and emis- 
sivity index. ISchmoll et al.l ()2009D apply the constraint 
that the inner emissivity index q\ should not exceed 5; if 
we impose the same constraint, we recover a very similar 
intermediate spin and inclination constraint. However, 
our base model allows the emissivity index to be much 
steeper (formally pegging at q\ = 10) and we find, as a 
result, that the inferred black hole spin tends to high 
values. For Suzaku B, our fits are in agreement with 
iPatrick et al. (2011) who finds that spin is unconstrained 
by this dataset. For XMM-Newton B, we find a spin pa- 
rameter and inc lination that is strongly discrep ant with 
those derived in lEmmanoulopoulos et all (|2011|) . We at- 
tribute this difference to the inclusion of a photoemission 
component, a newer absorption model (TBnew vs. wabs) 
and a different model construction. 

However, just looking at the results of our own uniform 
analysis, we find important and informative inconsisten- 

12 http://space.mit.edu/cxc/isis 
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Figure 5. Spectra, model [green] and residuals for pointing 
Suzaku-A (XIS0 [black], XIS1 [red], XIS3 [blue], PIN [purple]) for 
the fit with the base model. 



cies. Comparing the fit parameters derived from the dif- 
ferent observations we find that the black hole spin, the 
disk inclination, and the iron abundance appear to have 
significantly different values between the three datasets 
(Table [2]). Naturally we would expect these physical 
quantities to be constant on the time spanned by these 
observations. 

As we show below, a multi-epoch analysis reveals the 
need for a new continuum component to describe the 
soft excess. Foreshadowing that discussion, we note 
that adding such a soft excess component to the model 
for each individual dataset results in unconstrained spin 
parameters and emissivity indices, i.e., the individual 
datasets do not possess the S/N to determine spin and 
inclination in the more complex models; there is a signif- 
icant trade off between the parameters of the soft com- 
ponent and the spin and inclination for each individual 
data set. This also implies that that the extreme values of 
spin and q\ found in a fit of the base model to Suzaku-A 
may arise due to a model which is (possibly falsely) at- 
tempting to fit a very smooth soft excess with an ionized 
reflection spectrum, necessitating extreme broadening. 

3.3. Multi- epoch analysis 

The inconsistencies resulting from independently fit- 
ting our base model to the individual pointings leads 
us to the use of multi-epoch fitting. In multi-epoch fit- 
ting, we assume that the source spectrum is always com- 
posed of the same principal physical components and 
that, when physical considerations demand, parameters 
are forced to have the same value for all epochs. In 
our case, we demand that the black hole spin parameter, 
the inclination of the disk, the iron abundance, and the 
strength of the photoionizcd emission component (be- 
lieved to originate from an extended AGN wind) must 
have common values across the fits to all of the datasets. 
In all cases, these parameters are believed to remain 
fixed for any given AGN over human timescales. Aside 
from these additional constraints on the parameters, the 
model set-up is the same as for the fitting of the individ- 
ual pointings. We allow the slope/normalization of the 
primary powerlaw, as well as the normalization, emis- 
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sivity profile and ionization state of the disk reflection to 
float freely between datasets — these parameters depend 
upon the structure / geometry of the disk corona that can 
plausibly change on timescales that are (much) shorter 
than the inter-pointing spacing. Fitting in this manner 
we find that the broken power law emissivity profile is 
ill-constrained and a simple power law emissivity profile 
is sufficient to describe the data very well. Furthermore, 
to reduce the number of spectral bins in each dataset, 
we co-added the front-illuminated XIS data (XISO and 
XIS3); this eases the significant computational expenses 
of multi-epoch fitting. 

Statistically, fitting this model to the full multi-epoch 
dataset produces a good fit, with reduced \ 2 01 1-06. The 
best fit parameters and errors are listed in Table [3l The 
best fitting spin is a little higher (a — 0.7ll ; og ) and 
the inclination a little lower (i = 37j^ deg) than found 
in the previous study of ISchmoll et al.l (|2009f ) but, given 
the error bars, there is no strong discrepancy. The ion- 
ization parameter either takes rather low values (XMM, 
Suzaku A) or very high values (Suzaku B). This is di- 
rectly tied to the fact that, apart from the iron line 
itself, there are no other strong features in the soft X- 
ray band that the reflection models are locking onto 
in the spectrum. One aspect of this fit that gives us 
pause is the high PIN/XIS cross normalization compo- 
nent, Cxisf-PiN ~ 1-74. While we have already acknowl- 
edged that the pointing errors would render invalid the 
standard value of Cxisf-PiN = 1.18, this measured value 
would imply a ~ 30% suppression of the net XIS count 
rate below the fiducial value, a significantly larger de- 
pression than is reasonably expected from the pointing 
problems. A more likely possibility is that the spectrum 
has curvature in the concave sense so that the true PIN- 
band flux is higher than predicted by our base spectral 
model. 

Partly motivated by these considerations, but also 
guided by the known phenomenology of AGN, we explore 
the hypothesis that we are missing another component 
that can contribute to the soft excess (above and be- 
yond any soft excess emission from the blurred uniformly 
ionized reflection, that models the broad iron line). A 
common phenomenological model for the soft excess is a 
simple blackbody component. Adding a blackbody com- 
ponent to the spectral model does not yield improved 
fits — none of the datasets show the strong curvature 
representative of the Wien tail of a blackbody. How- 
ever, significant improvements in the goodness of fit are 
achieved using a soft excess model consisting of a "warm" 
thermal Comptonization component. In detail, we add 
a thermal Co mptonization com ponent described by the 
comptt model (Titarchuk 1994) with a seed photon tem- 
perature of 40 eV (representative of the expected thermal 
emission from the AGN disk), an optical depth r and an 
electron temperature kT. 

When applied to the multi-epoch data, the inclusion 
of this soft excess component leads to an improvement 
in the goodness of fit by A^ 2 = 198 for nine additional 
parameters. The best fit parameters and correspond- 
ing errors can be found in Table |4] As compared with 
the base model above, the best-fit spin has slightly de- 
creased to a — 0.52^q']J and the best- fit inclination has 
increased to i — 48j^ deg. While the change in spin 



between the two models is within the error bars, the in- 
clination change is significant and brings th e resul ts into 
line with the conclusions of ISchmoll et al.l ([2009). We 
find that the ionization state of the inner accretion disk 
is significantly higher at lower flux states. The biggest 
surprise is that, once we include the explicit soft ex- 
cess component, the iron abundance of the accretion disk 
is pushed towards the maximum value tabulated in the 
model, Z = IQZq. The formal 90% lower limit on iron 
abundance is Z > 8.2Z Q . From a modeling perspective, 
such a high iron abundance means simply that the iron 
line is strong in comparison to the Compton hump. 

However, another possibility is that the reflection from 
the inner accretion dis k is characterized b y multiple ion- 
ization components ()Nardini et al.l 1201 lh and a highly 
blurred, high-ionization reflection component can ac- 
count for the soft excess in a similar fashion than the 
additional Comptonization component. To explore this 
model, we refit the multi-epoch data with two disk re- 
flection components, each of which has its own ioniza- 
tion parameter and emissivity index. The only other 
difference to the multi-epoch fit above is that the iron 
abundance is assumed to be the same for all reflection 
components included in the model. Such a model leads 
to a fit solution of almost comparable quality than that 
employing the Comptonization-based soft excess compo- 
nent (Table [5]). In order to describe the smooth soft 
excess, one of the disk reflection components adopts a 
high ionization state (£ a ~ 1000) and a very steep emis- 
sivity profile (q a > 8) and the black hole spin parameter 
becomes large (a > 0.96). At the same time, much less 
broadening/blurring is needed in order to describe the 
broad iron line, and so the other disk reflection compo- 
nent adopts a low ionization (£& ~ 1 — 10) and a shallow 
emissivity profile (q < 2.2). This solution avoids the ex- 
treme iron abundance, with an inferred iron abundance 
of 0.68-0.85 Z Q . 

4. SUMMARY AND DISCUSSION 

The X-ray spectrum of Fairall 9 is one of the clean- 
est known in terms of having no discernible intrinsic ab- 
sorption. Fairall 9 has also been the subject of multi- 
ple pointings by XMM-Newton and Suzaku. This makes 
it an interesting target for a detailed study of the rela- 
tivistic reflection features including an assessment of the 
robustness of the derived parameters with respect to as- 
sumptions and choices made in the spectral modeling. 

A simple application of our base spectral model to the 
individual XMM-Newton and Suzaku pointings produces 
unphysical results, with inconsistent values for black hole 
spin, disk inclination and iron abundance being derived 
for each pointing. We identify several factors that are 
at play. We clearly see (and hence include in our base 
model) high ionization, narrow iron emission lines, pre- 
sumably from optically-thin photoionized plasma in the 
circumnuclear environment. When viewed with CCD- 
resolution (E/AE ~ 50 — 100), these narrow lines can 
blend with the blue wing of the broad iron line from 
the accretion disk, inducing strong degeneracies between 
the strength of the photoionized emission component and 
the accretion disk parameters. This degeneracy is par- 
ticularly strong for objects with intermediate inclination 
(35-60 deg) since, in these cases, the blue wing/edge of 
the iron line overlays the FeXXV Ka line complex. 
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Figure 6. Best fitting multi-epoch model (including the additional soft continuum source) for the Suzaku A (left), Suzaku B (middle) and 
XMM-Newton B (right) observations. Shown here is the total model (solid-thick black line), the primary power-law continuum (dashed 
black line), warm Comptonization component (dashed red line), the relativistically blurred ionized disk reflection (thick solid blue line), 
the distant reflector (thin solid magenta line), and the photoionized emission component (green line). 



To cut through this degeneracy, we make the assump- 
tion that the photoionized emission component is con- 
stant over the three years that span observations Suzaku- 
A, XMM-Newton-B, and Suzaku-B. We then employ a 
multi-epoch analysis, fitting our spectral model to these 
three datasets simultaneously — the black hole spin, ac- 
cretion disk inclination, iron abundance, and photoion- 
ized emission component are assumed constant across all 
three pointings, but all other spectral parameters and 
normalizations are allowed to float freely between point- 
ings. Since there is, indeed, significant flux variability 
between the three pointings, this approach is remarkably 
successful at breaking the degeneracy with the photoion- 
ized emission component and allowing us to constrain 
model values for the black hole spin and disk inclination. 
Of course, it must be highlighted that the constancy of 
the photoionized emission component is an assumption 
which could be violated if the emitting wind is compact 
(sub-parsec). This assumption must be tested by future 
high-resolution observations using the Chandra High En- 
ergy Transmission Gratings (HETG) or the Astro-H Soft 
X-ray Spectrometer (SXS). 

In addition, we find that the multi-epoch spectral fit is 
significantly improved by the addition of another spectral 
component. Our analysis therefore highlights the practi- 
cal importance of modeling the soft excess. The nature 
of the soft excess is still a matter of debate, and the soft 
Comptonization component that we attempt first is just 
one of the possible models discussed in the context of 
AGN spectra. In the case of Fairall 9, we have already 
excluded a soft excess that arises solely from inner disk 
reflection with a single ionization parameter, and have 
shown that an additional blackbody component is also 
rejected. 

We first consider the scenario in which the soft ex- 
cess is identified with an additional Comptonization 
component. As explicitly shown in Fig. [6j this addi- 
tional Comptonization component is the main contrib- 
utor to the soft excess during our two Suzaku observa- 
tions while ionized disk reflection dominates during the 
XMM-Newton B pointing. The inclusion of this soft 
continuum has a small effect on the best fitting black 
hole spin and disk inclination (with our final values be- 
ing a = 0.52^ 15 an d i = 48 deg), but has a dramatic 



effect on the inferred iron abundance of the accretion 
disk (Z > 8. 2Z(7)). With an expected abundance va lue of 
only 2-3Z Q ([G roves. Heckman fc Kauff mann 2006), such 
a high iron abundance implies an unusual star formation 
history in the galactic core, or some process that pref- 
erentially differentiates iron into the photosphere of the 
inner disk (e.g., see discussion of Reynolds et al. 2012). 

It is interesting that, in our fiducial multi-epoch fit 
(Table U, the Suzaku-B dataset requires a high emissiv- 
ity index, q > 6. Such high emissivity indices are nor- 
mally attributed to the action of extreme light bending 
focusing X-rays onto the innermost regions of the disk 
or, alternatively, the dissipation of work done by strong 
torques at the ISCO. Both of these explanations would 
require a black hole that is spinning rapidly (a > 0.9), as 
opposed to the inte rmediate spin (a « 0.5) that we infer 
(jFabian et al .1 1201 2). A possible resolution of this contra- 
diction comes from the vertical geometry of the accretion 
disk. Fairall 9 is a lumino us source accreting w ith an Ed- 
dington rati o of £ ~ 0.15 (lLohfink et al. 2012). Standard 
disk theory IShakura fc Sunvaevl (| 19731) tells us that the 
inner regions of the accretion disk will be very radiation 
pressure dominated and, away from the inner bound- 
ary, will have a disk thickness h ~ (3£/2?7)r 9 w 2.3r g 
where we have taken the radiative efficiency of the disk 
to be r] = 0.1. As the disk approaches the ISCO (at 
Hsco ~ 4.5r ff for a — 0.5) the thickness diminishes and, 
within the ISCO, the accretion flow forms a thin sheet 
that spirals into the black hole. This gives the photo- 
sphere of the inner disk the geometry of shallow bowl. If 
the X-ray source is located very close to the black hole 
(e.g. in a spin-powered magnetosphere) or on the inward 
facing surface of this bowl, the region close to the ISCO 
can be strongly irradiated while the disk surface at large 
distance may receive very little irradiation — this would 
manifest itself in our spectral analysis as a steep emissiv- 
ity profile. Alternatively, the high emissivity index may 
be a mirage induced by a steep iron abundance gradi- 
ent. Specifically, if the high iron abundance is due to 
radiative-levitation in the disk photosphere (?), it will 
have a strong radial gradient which would appear in our 
fits as a steep emissivity index. 

Alternatively, this steep emissivity profile could be the 
artifact of incorrect modeling of contamination in the 
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Figure 7. A comparison of the two best fitting spectral models, the double reflection model (blue) and the one with additional Comp- 
tonization soft excess (red), for the first Suzaku pointing (left) and the XMM pointing (right). The models differ significantly in the Nustar 
energy range marked with the black bar. 



spectrum. As outlined in § 12.1. 11 this Suzaku pointing is 
certainly affected by strong contamination in XISO below 
1 keV. We assumed for this analysis that the other two 
detectors were correct below 1 keV, but the steep pro- 
file could be an indication that this assumption is not 
justified. 

The alternative scenario that we consider for the soft 
excess in Fairall 9 is a model consistent of two ionized 
reflection components with different ionization states. 
This double-ionized reflection solution does present dif- 
ficulties, however, when we attempt to interpret it with 
physical models. There are two classes of model in which 
multiple ionization states may be present in the disk re- 
flection spectrum. Firstly, at any given radius, the den- 
sity of the disk photosphere may be patchy leading to 
mixed regions of high and low ionization. In such a sce- 
nario, we would expect the two ionization components 
to share approximately the same emissivity index since 
they respond to the same irradiation profile. This dis- 
agrees with the fact that we see very different emissiv- 
ity indices between these two components and allows us 
to reject this patchy disk model. Secondly, we may ex- 
pect the ionization parameter of the disk surface to have 
a strong radial gradient due to the radial dependence 
of both the irradiating flux and the disk density. The 
fact that we see a high-ionization component which is 
highly blurred (coming from the smallest radii) and a 
low-ionization component which is more weakly blurred 
(coming from substantially larger radii) fits in with this 
picture. However, preliminary theoretical calculations of 
the reflection from a disk with a radial ionization gradi- 
ent reveals that the intermediate ionization states should 
also imprint strong soft X-ray features which we do not 



see (|Garcia et al.H2012D . 

Fundamentally, however, it is difficult to distinguish 
these two scenarios (new Comptonization continuum or 
second disk reflection component) with existing data. Al- 
though the models do not differ significantly in quality 
of their description of the data, the physical interpreta- 
tions are quite distinct. Plotting the models to energies 
well above the PIN energy range (Fig. [7]) points out how 
future hard X-ray studies with Nustar or Astro- H could 
break this degeneracy. The hard X-ray flux predicted 
by the double reflection scenario is substantially smaller 
than the one predicted by the model including the addi- 
tional Comptonization component. In addition, the high 
spectral resolution will help in the soft X-ray part putting 
firm constraints on the smoothness of the X-ray excess. 

Fairall 9, with its absence of warm absorption, is an in- 
teresting target to study the properties of the soft excess. 
This enables us to test possible scenarios explaining the 
soft excess in the spectrum, as well as the time variability 
of this component. We note that the shape of the addi- 
tional soft excess seems to change as the source becomes 
dimmer (Fig [5]). The dimmer the source the steeper the 
additional component. At the time the contribution of 
the additional component decreases and the lower ion- 
ized reflection dominates. A more detailed analysis of 
this phenomenology could lead to more insight into the 
nature of this additional component but requires more 
high quality pointings and is therefore subject to future 
studies. 

As a final comment, we highlight that the Suzaku-H 
pointing needs to treated carefully and the cross nor- 
malization between XIS and HXD-PIN cannot be ex- 
pected to be at normal value but instead is found to 
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be higher (1.3). We n ote that the reported discovery 
by iTatum et al.l ([2012) of a new spectral component 
that only reveals itself in the PIN-band of this partic- 
ular dataset (and interpreted as a transitory marginally 
Compton-thick absorber) may be affected by this incor- 
rect cross-normalization. 
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Table 2 

Spectral Parameters for Fairall 9 base model fits with additional photoionized emission; See the text for a description of the model. 
Suzaku Spectra are normalized to XISO data. The power law normalization is photons kcV -1 cm~ 2 s — 1 at IkeV. 
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Table 3 

Spectral Parameters for Fairall 9 multi-epoch base model fit with additional photoionized emission; See the text for a description of the 
model. Suzaku Spectra are normalized to XISO data. The power law normalization is photons keV -1 cm~ 2 s _1 at IkeV. 
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Table 4 

Spectral Parameters for Fairall 9 multi-epoch fit with additional photoionized emission and soft excess component; See the text for a 
description of the model. Suzaku Spectra are normalized to XISO data. The power law normalization is photons keV - 1 cm~ 2 s -1 at 1 keV. 







Suzaku A 


Suzaku B 


XMM B 


contiuum & 
cold reflection 


A pex [io- 3 ] 
r 

R 


6 71+ u + 3 
D - ' -"--0.78 

1 90+ 04 

i- yu _o.os 
0.84+°'°! 

—0.07 


3.59ti;» 

i 7-.+0.15 
1 -' 1 -0.09 

0.79+°" 

—0.17 


3 60+"' 17 
J - DU _0.18 

1 85+ ' 02 
l- oo -0.02 

1.08+°'" 

—0.11 


ionized reflection 


Arfionx [IO" 4 ] 

Fe/ Solar 


0.57+^ 


99+ ua4 
10+? 77 

— 1.77 


13+0-02 




£ [erg cms -1 ] 


20.9+^ 


lO.lt^ 4 


or,, r- + 381.0 

3513.5I 393 7 


relativistic blurring 


1 

i [uegj 


o +0.23 
Z -0 


9.5t°i 
52+ ' 19 
48+ 6 


O+0.33 


Comptonization 


A M n— 31 
^-compTT [10 J 

HI [KCVJ 

T 


A m +u.yu 

4 - 01 -1.88 
on ko+49.09 
20.52_ 953 

0.52l°;g 


j- 17+U.Ul 

°- i '-0.95 
or ns+ 1 - 40 

~ 07 


q « K +18b2b.b3 
3 - 65 -1.65 

no+11.29 
°°- zo -33 .28 
0.01+Q 


plasma 


Aphot 

log? 




6 OO+0.67 




cross calibration 


CxiSf 
CXIS1 
c XISf-PIN 


1.00 
1.00 
1.16 


1.00 

oo+o-oi 

1 26+ ' 07 




individual 


X 2 /dof 


2232.56/[2106-ll] 


1985.06/[2060-13] 


1128.64/[1100-11] 


contributions 


p- value 


0.02 


0.83 


0.20 




X 2 /dof 




5346.275/ [5276-34] 






p- value 




0.16 










1.02 





Table 5 

Spectral Parameters for Fairall 9 multi-epoch fit with additional photoionized emission and two ionized reflection components; Sec the 
text for a description of the model. Suzaku Spectra are normalized to XISO data. The power law normalization is 

photons keV" 1 cm~ 2 s -i at IkeV. 







Suzaku A 


Suzaku B 


XMM B 


contiuum & 
cold reflection 


A pcx [IO" 3 ] 

r 

R 


8.02t U :& 
2.07+.°;°! 

1 OO + 0.21 

1 - JO -0.19 


9.84±°'j< 

9 , 7 +0.01 
z - ± ' -0.01 
n qc+0.23 
U.95_ 16 


4.81+1^ 

2 09+ ' 04 
z - ua -0.04 

1 89+ ' 11 
± - oa -0.16 


ionized reflection 1 


A rcflionx [IO" 4 ] 


51+0-00 
u - oi -0.12 


u - 41 -0.28 


U -0 




Fe/ Solar 
£ [ergems - *] 




0.751°+° 
1091.9+™ 


100 +6900 


relativistic blurring 1 


q 

a 

i [deg] 


io.oIH 


io.o±r 3 

n Q-+0.02 
"•"'-0.01 

36+ 3 


10.0™ 


ionized reflection 2 


A rcflionx [io- 4 ] 

£ [crgcms~i] 


1 00+°'^ 
!- uu -0.35 
19 1+8-6 
1^-1-5.9 


o or+o.w 

Z.OJ_ 43 
-1+1.2 
1-0 


i CT9+0.39 
!- JZ -1.02 

1+30.6 

1-0 


relativistic blurring 2 


i 


9+0.2 

z -0 


9 i+U.3 
z -l-0.1 


5.6±+ 4 


plasma 


Aphot 

log? 




0-09+°; 3 ° 

4 O+0.7 




cross calibration 


CxiSf 
CXIS1 
c XISf-PIN 


1.00 
1.00 
1.16 


1.00 

gg+O-Ol 

1 43+ - 07 




individual 


X 2 /dof 


2208.44/[2106-ll] 


2012.66/[2060-13] 


1168.47/[1100-11] 


contributions 


p- value 


0.04 


0.67 


0.05 




x 2 /dof 




5389.6/ [5276-34] 






p- value 




0.08 










1.02 





